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Abstract 



We calculate one-loop quantum energies in a renormalizable self- 
interacting theory in one spatial dimension by summing the zero-point en- 
ergies of small oscillations around a classical field configuration, which need 
not be a solution of the classical field equations. We unambiguously imple- 
ment standard perturbative renormalization using phase shifts and the Born 
approximation. We illustrate our method by calculating the quantum energy 
of a soliton/antisoliton pair as a function of their separation. This energy 
includes an imaginary part that gives a quantum decay rate and is associ- 
ated with a level crossing in the solutions to the classical field equation in the 
presence of the source that maintains the soliton/antisoliton pair. 
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I. INTRODUCTION 



The study of exactly soluble 1+1 dimensional problems has yielded many insights into 
fundamental problems in field theory. Other 1+1 dimensional problems cannot be solved 
exactly, making it important to understand which properties of exact results will generalize 
to more generic cases and which are special to the exactly soluble cases. A particularly 
rich area is the study of solitons — time independent non-linear solutions — that arise is 
interesting theories in both one and three dimensions. Renormalized quantum fluctuations 
can modify a soliton's properties, perhaps even destabilizing a classically stable solution 
or stabilizing a classically unstable one. The renormalization process is simpler in one 
dimension than it is in three, but it still must be approached carefully. Merely "canceling 
the infinities" is clearly not sufficient when we want to regard a finite result as a physical 
prediction of a particular theory defined under fixed renormalization conditions. 

In this Letter we speciahze the discussion of [1] to include the case of 1-dimensional 
self-interacting theories. We develop an unambiguous procedure for applying standard per- 
turbative renormalization methods to the calculation of the one-loop "effective energy" of a 
time-independent, spatially symmetric field configuration. Our method is also practical and 
efficient for numerical calculation, which we illustrate by considering a time-independent 
family of configurations (f){x, Xq) that continuously interpolates between the trivial vacuum 
{xq — 0) and a soliton/antisoliton pair (xq — > oo) in 0^ theory. The energy of these con- 
figiuations can include an imaginary part, which we show measures the instability of the 
system to quantum fiuctuations. We also consider the exactly soluble cases such as the 
xo ^ oo limit of the soliton/antisoliton pair, and obtain standard results that are free of 
cutoff ambiguities. Throughout, we will be able to gain insight into these field-theoretic 
problems by understanding them in terms of physical quantities like phase shifts and bound 
state energies. 

The present work can be regarded in several ways. On the one hand, one-dimensional 
models provide very clean pedagogical examples of the general methods developed in [1]. On 
the other hand, quantum energies of one dimensional solitons are interesting in their own 
right, as [2], [3], and [4] attest. Our methods allow us to resolve questions that are obscure 
if one can only analyze quantum fluctuations about solutions to the classical equations of 
motion. Our results on the quantum fluctuations about the 0^ kink are not new — they 
coincide with the well known results of Dashcn, Hasslachcr and Ncvcu [5]. However, our 
method displays interesting features like the correlation between imaginary frequency small 
oscillations and level crossing in the space of solutions to the equations of motion, and 
provides a simple method for avoiding ambiguities. 



To illustrate our method, we will consider a standard 0"^ theory with spontaneous sym- 
metry breaking and a source J[x). The action is 



II. FORMALISM 




(1) 



2 



where is the countcrtcrm Lagrangian. We have rescaled the field somewhat uncon- 
ventionally in order to make explicit the correspondence between the powers of A and the 
powers of U. Classical terms will go as the one-loop terms we compute will go as A°, and 
higher loops will contribute with higher powers of A. The mass of fiuctuations around the 
trivial vacua (t){x) — ±1 is m, and we define the potential U{(f)) — ^(<^^ — 1)^. 

We consider a fixed field configuration (l)o{x). Initially we will assume that (/»o(a;) = 
0o(— ^)) which restricts us to the topologically trivial sector of the theory. However, we will 
see that our method works equally well for the case of (pQ^x) = —4>o{—x), so that in fact we 
can deal with configurations with any topology as long as ?7(0o) has reflection symmetry. 
We adjust the source so that (j)o is a stationary point of the action, which means that J then 
solves the equation 

J{x)^^-lm^{<Pl-cf>o). (2) 

We can always solve this equation for J, but there is no guarantee that the J we find will 
always correspond to a unique ^o, as we will see later. If is a solution to the equations of 
motion, of course the source will be zero. 

We would like to consider the leading quantum correction to the classical energy of this 
configuration, which we can represent as the sum of the zero-point energies of the normal 
modes of small oscillations around 0o- Writing (p = (pQ + rj, the normal modes are solutions 
of 

-^ + {V{x) + m^)ri^co'^r] (3) 

with V{x) = U"{4>q{x)) — rn? = |m^(0o(^) ~ !)• The quantum change in energy in going 
from the trivial vacuum to ^o{x) is then 

^[M ^\{Y.uj-Y.^o)+S^^^E + £^ (4) 

where ujq are the free solutions (with (j)'^{x) — 1), and S^t is the contribution from the 
counterterms. 

Note that it is possible that some of the values of uP' will be negative, so that in these 
directions our stationary point is a local maximum rather than a minimum of the action. 
These solutions will add an imaginary part to the energy, which we can interpret via analytic 
continuation as giving the decay rates through the unstable modes [6]. If there is a direction 
in field space in which small oscillations lower the energy, we should be able to keep going 
in that direction and arrive at a lower minimum. Later, we will see explicitly that the 
appearance of unstable modes coincides with the existence of a second solution to the 0o 
equation with the same J and lower energy. 

We would like to rewrite the sum over zero-point energies as an integral over phase space 
of the product of the energy and the density of states. We can then break this integral 
into a sum over bound states and an integral over a continuum, representing the latter 
in terms of phase shifts. In order to do so, however, we must review the peculiarities of 
Levinson's theorem in one dimension. For more details on these results see [7] and [8]. For 
a symmetric V{x), we can divide the continuum states into symmetric and antisymmetric 
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channels, and then calculate the phase shift as a function of k separately for each channel, 
where u"^ = k"^ + m^. The antisymmetric channel is completely equivalent to the I = case 
in three dimensions, so we have 

5a(0) = nATT (5) 

where is the number of antisymmetric bound states. However, we must be careful in 
dealing with the special case of a state exactly at = 0. In this case the solution to eq. (3) 
with A: = goes asymptotically to a constant clS X ^ OO ^ clS opposed to the generic case 
where the /c = solution goes to a constant plus linear terms in x. Just as in the / = 
case in three dimensions, this state contributes \ to n^- We will refer to such states as 
"half-bound states." 

In the symmetric channel, Levinson's theorem becomes 

5s(0)=ns7r-| (6) 

where a bound state at A; = contributes | to ng, just as in the antisymmetric case. We 
can see the importance of getting the half-bound states right by looking at the free case: 
then the phase shift is zero everywhere, and the right-hand side is zero because the free case 
has a half-bound state (the wavefunction = constant). The situation is equally subtle for 
reflectionless potentials, all of which have half-bound states and 5s (0) + 5a(0) equal to an 
integer times tt. 

We are now ready to rewrite the change in the zero-point energies in terms of phase 
shifts. Letting Ej be the bound state energies (again with /c = bound states contributing 
with a I), p{k) be the density of states and po{k) be the free density of states, we have 



where uj{k) — y/W+lr^ and we have used 

p{k)=Poik) + ^^{SA{k) + Ss{k)). (8) 

Note that the ^ term subtracts the contribution of the free half-bound state. 

Eq. (7) is divergent (the phase shifts fall as | for A; — > oo), which is what we should expect 
since it includes the divergent contribution from the tadpole graph without the divergent 
contribution from the counterterms that cancels it. 

To avoid infrared problems later, we first use Levinson's theorem to compute the change 
in particle number, which is given by 







where again, half-bound states are counted with a | in the sum over j and the — | comes 
from the contribution of the free half-bound state. Subtracting y times this equation from 
eq. (7), we have 
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^E^- Y.{E, -m) + l -{u{k) - m)-{6^{k) + Ss{k)). (10) 

Next, we follow [1] and subtract the first Born approximation to eq. (10), which corre- 
sponds exactly to the tadpole graph. We must then add it back in using ordinary renormal- 

ized perturbation theory. However, in 1+1 dimensions we can adopt the simple renormal- 
ization condition that the countertcrms cancel the tadpole graph and perform no additional 
finite renormalizations beyond this cancellation. With this choice, there is then nothing to 
add back in. 

The first Born approximation is given by 

(k) = —— / V (x) cos^ kx dx 

K Jo 

S^^\k) = --j^ V{x)sm^ kxdx. (11) 

Notice that the sum of these two depends on V{x) only through the quantity {V) — 
dxV{x), so we can indeed cancel the tadpole contribution with available countertcrms. 
Subtracting the first Born approximation, we have 

1 r°° dk d iV\ 

^[M = 2 E(^. - H + i, ^(^(^) - ^)^(^^(^) + ^s(^) + ^^^^ 

which is completely finite and well-defined, since the integrand goes like p- for — > oo and 
goes to a constant at A: = 0. In particular, we are free to integrate by parts, giving an 
expression that will be easier to deal with computationally, 

1 dk k (V) 

£M--,UE,-m)-l + W + V). (13) 

The subtraction has canceled the leading | behavior of the phase shift, leaving an integral 
that is finite in the ultraviolet. It must, since in 1+1 dimensions, eliminating the tadpole is 
sufficient to render the theory finite. (In 3+1 dimensions, we can decompose the problem 
into a sum of 1-dimcnsional problems, labeled by the angular momentum £. We then need 
two subtractions: the first subtraction renders each 1-dimensional problem finite, and the 
second renders the sum over i finite.) 

Our use of Levinson's theorem to "subtract" AE in the infrared (replacing eq. (7) by 
eq. (10)) avoided a subtlety of the Born approximation in one dimension that does not 
occur in higher dimensions: as A; — > 0, the symmetric contribution from eq. (11) introduces 
a spurious infrared divergence, since it goes like i. (Each higher dimension adds a power of 
k near A; = 0. We can see this in the Feynman diagram calculation, where these powers of 
k come from the measure.) 

III. APPLICATIONS 

We can now use eq. (13) to calculate quantum energies for specific field configurations. 
Knowing that our model has a "kink" soliton solution (f)Q{x) — tanh^, we will consider a 
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family of field configurations that continuously interpolate between the trivial vacuum and 
a soliton/antisoliton pair, 

4>q{x, xq) — tanh —{x + xq) — tanh —{x — xq) — 1, (14) 

with 2xq measuring approximately the separation between the soliton and antisoliton. Unlike 
the kink, these configurations are not solutions of the equations of motion, except in the 
Xq ^ Q and Xq ^ oo limits. Thus we will need to to introduce a source that will vanish in 
these limits, and we must analyze the stability questions we raised earlier. 

For xq very small, we simply have a small attractive perturbation from the trivial vacuum 
held in place by a small source, which we would not expect to introduce any instabilities. 
In terms of the scattering problem, for small Xq, the potential is too weak to bind a state 
with a binding energy greater than m, which would give an imaginary eigenvalue. 

For xq very large, we have a widely separated soliton/antisoliton pair. We know from 
translation invariance that a single soliton has a mode with uP' — Since this zero mode 
corresponds to a nodclcss wavefunction, it is the lowest energy mode. The soliton/antisoliton 
pair has two translation modes that will mix, giving a symmetric eigcnstate with a slightly 
lower energy and an antisymmetric eigenstate with a slightly higher energy. Thus we expect 
to find a single symmetric mode with oJ^ < 0, which will contribute an imaginary part to 
£^[(?i>o]. The imaginary part gives the rate for our field configuration to decay through this 
mode toward the trivial vacuum. As xq — oo, we should find that the imaginary part goes 
to zero, and the real part goes to twice the energy of a single soliton, which we compute 
exactly below. 

According to this analysis, there should be a finite, nonzero value of xq where the imag- 
inary eigenvalue first appears. At this point, the field becomes unstable with respect to 
small perturbations in some direction in field space. Therefore, the energy must have a 
lower minimum that we can reach by moving in that direction. This configuration ^fJo{x, xq) 
is a second stationary point of the action with the same source, which crosses 0o at the 
value of Xq where the imaginary part for the energy appears. (For smaller values of Xq, this 
solution still exists but has higher energy.) In terms of the scattering problem, this crossing 
appears when the potential has a bound state with u'^ = 0. We can identify this state 
explicitly: since (po and i/jq satisfy eq. (2) with the same J, as we approach the crossing, the 
wavefunction r]{x) = i/Jo{X; x^) — <po{x, Xq) becomes a solution to eq. (3) with uj"^ = 0. 

We have carried out this computation numerically, and find results that agree with all 
of these expectations. To compute the antisymmetric phase shift, we use the method of [1] 
to write 

S^ik) = -2ReP{k,0) (15) 
where P{k, x) is at a; = oo and satisfies 

+ 2kp' + 2{f5f + V(x) = 0, (16) 

with prime denoting differentiation with respect to x. To obtain the symmetric phase shift, 
we can follow the same derivation, but instead of imposing -0(0) = on the wavefunction, 
we instead impose V''(0) = 0, giving 
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e^'^s = ' ^ (17) 

e-^'f^* k + 2(3'* ^ ' 



with j3 the same as above, so that 



^s(fc) = ^AW-arctan(^^^^^i^). (18) 
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FIG. 1. The real and imaginary parts of <?[(/'o], with £ in units of m, for xq increasing from 
to ^. 

m 

Figure 1 shows the real and imaginary parts of £[(t)()\ as a functions of Xq, starting from 
the origin at xq = 0. When xq ~ a single imaginary eigenvalue appears and S[(f)o] 
becomes complex. For xq large, ^[0o] approaches 2(^^ — ^), twice the standard result for 
a single soliton. Figure 2 shows the path of the energy in the complex plane parametrically 
as a function of the separation. The actual trajectory of S[(j)Q] in the complex plane has 
little significance, since it depends in detail on the functional form of (f){x, Xq). However, the 
general features — beginning at the origin, moving up the real axis, out into the complex 
plane, and finally asymptotically to the real two-solution value — are characteristics of any 
00 that begins at the vacuum and ends at a well separated kink and antikink. 

Fig. 3 shows (po{x, xq) and the second solution to the equations of motion with the same 
J, i/jq^XjXo). ipo goes to the trivial vacuum as xq — > oo, and becomes a widely separated 
soliton/antisoliton pair as xq 0, crossing 0o at xq ~ 0.75. For xq below this value, 0o has 
lower classical energy, while above this point, i/jq has lower classical energy, and this crossing 
appears precisely where the imaginary eigenvalue appears in the small oscillations spectrum. 

Through this continuous deformation from the trivial vacuum, we have arrived at a 
widely separated soliton/antisoliton pair, which we can now separate into independent con- 
figurations with nontrivial topology. These configurations are exactly soluble, so we will be 
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FIG. 2. The trajectory of £[(j)o] in the complex plane, with £ in units of m, for xq increasing 
from to in steps of 0.0625/m. 



able to study them analytically. However, our method does not rely on having an analytic 
solution, so we could numerically calculate the energy of a generic field configuration with 
nontrivial topology using the same techniques. 

In one dimensional quantum mechanics, potentials of the form 



VAx) 



"h 1 2 1 9 / ^^-^ \ 

— — m seen {——) 



(19) 



with £ an integer are exactly soluble and refiectionless. The single soliton solution in our 
model, <j)o{x) = tanh corresponds to V{x) = — |m^sech^^, the i = 2 case of this family. 
(The sine-Gordon soliton corresponds to the £ = 1 case.) For a refiectionless potential, 
Ss{k) = 5a (A;); to reconcile this equality with Levinson's theorem in the symmetric and 
antisymmetric channels, Vi{x) must have a half-bound state. (We saw this behavior already 
for the i = case, the free particle.) We also note that although 6s = 5a, ^g^'' 7^ 5a\ so 
the renormalized contributions from the symmetric and antisymmetric channels are not the 
same. In addition, the bound state contributions will also be unequal. 
For our case, the exact result for the phase shift is 



Ss{k) = 5a(^) = — arctan 



3mk 



(20) 



where the branch of the arctangent is chosen so that the phase shift is continuous and goes 
to zero for k ^ 00. The Born approximation is 



6l^\k)+6^\k) 



3m 



There are three bound states: a translation mode with E = 0, a state with E 
and a half-bound state with E = m. Using eq. (12) or eq. (13) we find 



(21) 
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FIG. 3. Solutions to the equations of motion with a source given by eq. (2), as functions of x 
and xq in units of ^. The shaded graph is (j)o{x,xo), which is guaranteed to be a solution by the 
construction of J. The unshaded graph gives the second solution iIjq{x,xo). (A finer mesh is used 
for this graph between xq = and xq = ^ in order to illustrate its behavior in this region.) 



in agreement with [5], [4], and the mode number cutoff method in [2]. 

We can calculate the energy of the soliton of the sine-Gordon model using the same 
methods. In this model, the 0^ potential is replaced by U{(f)) = m'^{cos(p — 1), which has 
soliton (anti-soliton) solutions 0o = 4 arctan(exp(±m(x — xq)))- The phase shift is 

771 

Ss{k) = 5a(^) = arctan — (23) 

with Born approximation 

The bound states are just the translation mode at = and the half-bound state ai E = tu, 
giving 

Ttl 

£[<t>o] = -- (25) 

TT 

again agreeing with the established results. 
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IV. CONCLUSIONS 



We have demonstrated that our method allows us to calculate one-loop energies of a 
wide variety of static, symmetric field configurations. We have avoided ever having to deal 
directly with subtleties of the cutoff or the boundary conditions, relying instead on the 
formalism of the Born approximation and Levinson's theorem. Using this formahsm, we can 
apply renormahzation conditions fixed in the topologically trivial sector so that our theory 
gives unambiguous predictions in terms of well-defined masses and coupling constants, and 
express the result in terms of simple physical quantities that allow for robust numerical and 
analytic calculation. 
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